NON-UNIFORM BOUNDS IN LOCAL LIMIT THEOREMS 
IN CASE OF FRACTIONAL MOMENTS 



S. G. BOBKOV^^*, G. P. CHISTYAKOV^-'', AND F. GOTZE'*''* 

Abstract. Edgeworth-type expansions for convolutions of probability densities and powers 
of the characteristic functions with non-uniform error terms are established for i.i.d. random 
variables with finite (fractional) moments of order s > 2, where s may be noninteger. 



1. Introduction 

Let {Xn)n>i be independent and identically distributed random viariables with 'EXi = 
and BXf = l7 

If Xi has finite moments of all orders, and if the densities pn of the normalized sums 
Sn = {Xi + • • • + Xn)/\/n exist, they admit a formal Edgeworth-type expansion with respect 
to the powers of l/\/ra 

oo 

p„(x) = ^(x)+^gfc(x)n-*^/2. (1.1) 

k=l 

Here, ip{x) = e~^^/^ denotes the density of the standard normal law. 



*w = ,wE'^.«.Wj;^(|)"...(^)". (1.2) 

where Hk{x) are the Chebyshev-Hermite polynomials, and 

ik 

denote the cumulants of the underlying distribution. The summation in (1.2) runs over all 
non-negative integer solutions (pi, . . . -.Pk) to the equation pi + 2p2 + • • • + kp^ = k with 
j =Pl-\ hPfc- 

A precise asymptotic statement about the formal series (1.1) requires that some moment 
E|Xi|'^ of order s > 2 is finite (while the moments of higher orders may be infinite). In this 
case, the k-th order cumulants are well-defined for the values k = 1, . . . ,m, and respectively, 
the functions qk are defined for k < ni — 2, where m = [s] is the integer part of s. Therefore, 
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one needs to evaluate the error of the approximation of pn by the fohowing partial sums of the 
series (1.1), 

m-2 

(Pm{x) = ip{x) + ^ qk{x) n"''/^, m = [s]. 
k=l 

One of the aims of this paper is to prove the following: 

Theorem 1.1. Assume that E|Xi|* < +oo, for some s >2. Suppose has a bounded 
density for some no ■ Then for all n large enough, Sn have continuous densities pn satis- 
fying, as n oo, 

(1 + IxD {pn{x) - ifmix)) = o(n-(^-2)/2) (1.3) 

uniformly for all x. Moreover, 
(1 + \xY){pn{x) - ^rn{x)) = o(n-(^-2)/2) + (l + |x|^-"*) (©(n-^^^-^V^) + o(n-(^-2))). (1.4) 

In fact, the implied sequence and constants in the error terms hold uniformly in the class 
of all densities with precribed moment tail function t — >■ E \X\Y l{|Xi|>t}5 parameter no and a 
bound on the density pno ■ 

For \x\ of order 1, or when s = m is integer, both relations are equivalent. But for large 
values of |x| and s > m, the assertion (1.4) gives an improvement over (1.3), which is essential 
in some applications. 

If 2 < s < 3, (1.4) becomes 

(1 + \x\'){pn{x) - ip{x)) = o(n-(^-2)/2^ ^ ^ |.x|^-2)o(n-(^-2)). 

In particular, for the smallest value s = 2, this contains the Gnedenko local limit theorem 
sup^, \pn{x) — ip{x)\ — 0, as n — oo. 

If s = m is integer and m > 3, Theorem 1.1 is well-known; (1.3)-(1.4) then simplify to 

(1 + \xn{Pnix) - Vm{x)) = o(n-("^-2)/2). (1.5) 

In this formulation the result is due to Petrov [Pel] (cf. also Pctrov [Pe2], p. 211, or Bhat- 
tacharya and Ranga Rao [B-RR], p. 192). Without the term 1 + \x\"'', (1.5) can be found in the 
classical book [G-K] ; this weaker variant goes back to the results by Gnedenko [G] and an ear- 
lier work by Cramer (who used, according to [G-K] , additional assumptions on the underlying 
density) . 

Thus, Theorem 1.1 extends these well-known results to the case, where s is not necessarily 
integer. The range 2 < s < 3 is of interest, as well. Our interest in these somewhat tech- 
nical extensions, especially (1.4), was motivated by open questions as to the actual rate of 
convergence in the so-called entropic central limit theorem. Here the relation (1.4) led to an 
unexpected behaviour of the error in the approximation of the entropy of sums of independent 
summands when s increases from 2 to 4. (This error stabilizes at s = 4, in contrast with the 
usual Berry-Esseen-type theorem for distribution functions, where stabilization of errors starts 
at s = 3). In the entropic central limit theorem the classical non- uniform bound (1.5) is not 
precise enough to derive upper bounds for errors. 
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Note that the assumption of boundedness of p„ in Theorem 1.1 (for some n or, equivalently, 
for all large n) is necessary for conclusions such as (1.3)-(1.5). It is equivalent to the property 
that the characteristic function v{t) = Ee**"^i is integrable with some power v > 1, i.e., 

/ + 00 
\v {t)\'' dt < +00. (1.6) 
-oo 

In this case pn are bounded for all n > 21/. The condition (1.6) is sometimes called "smooth- 
ness" ; it appears naturally in many problems of the asymptotic behaviour of the densities (see 
e.g. [S] for detailed discussion). 

Nevertheless, this condition may be removed at all, if we require that (1.3)-(1.4) hold true 
for slightly modified densities, rather than for p„. 

Theorem 1.2. Let E iXil'* < +oo, for some s > 2. Let c denote an arbitrary number with 
< c < 1 . Suppose that for n large enough Sn have absolutely continuous distributions with 
densities pn- Then, for some probability densities pn, 

a) The relations (1.3)-(1.4) hold true for Pn in place of Pn ; 

b) /J*^ \Pn{x) — Pn{x)\ dx < c^, for all n large enough; 

c) Pn{x) = Pn{x) almost everywhere, if Pn is bounded (a.e.) 

It seems that Theorem 1.2 has not been stated in the literature, even when s is integer. 
Here, the property c) is added to include Theorem 1.1 in Theorem 1.2 as a particular case. 

It turns out that the statement of Theorem 1.2 is more appropriate for a number of appli- 
cations. For example, it implies that pn — (fm ^ in the mean, i.e., there is convergence in 
total variation norm for the corresponding distributions with rate 

/+00 
\pn{x) - <^„(x)| dx = o(n-(^-2)/2). (1.7) 
-oo 

For s = 2 and 952(2:) = (^(x), this statement corresponds to a theorem of Prokhorov [Pr], while 
for s = 3 and (fsix) = (p{x)[l + as - to the result of Sirazhdinov and Mamatov [S-M] 

(they also covered the case 2 < s < 3 with O in place of in (1.7)). If s > 3 is integer, (1.7) is 
mentioned in [Pe2] for a more general L^-convergence — however, under the assumption that 
the densities p„ are bounded. 

Theorem 1.2 allows to study non-uniform convergence in (1.3)-(1.4) as well when excluding 
exceptional "small" sets (via additional assumptions of entropic type). 

The proof of Theorems 1.1 and 1.2, which are formally based on the application of the in- 
verse Fourier transforms, involves operators, namely, Liouville fractional integrals and deriva- 
tives. For this step, we analyse the decay of the Fourier transform for special classes of finite 
measures with finite fractional moments. (Apparently standard truncation methods are much 
to density-sensitive and do not provide the required asymptotics.) An essential part of the 
argument is devoted to the routine analysis of powers of the characteristic functions and more 
general Fourier transforms in Edgeworth-type expansions. For this step the requirement (1.6) 
is irrelevant. 
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In order to describe one of the main intermediate results, which is, as we beheve, of an 
independent interest, let us start with a random variable X, such that EX = 0, EX^ = 1, 
and E < +00, for some s > 2. Introduce the characteristic function v{t) = Ee**'''', t G R. 

If m = [s], the normalized powers Vn{t) = that is, the characteristic functions of 

Sn, can be approximated by the functions 



(m—2 N 
k=l ' 



Here we use the classical polynomials 



E 



1 



1L\ 



(fc + 2)! 



Vh 



^fc+2(pi+-+pfc) 



pi+2p2H |-A;pfc=A; ' 

of degree S/c, where the summation is performed as in (1.2). Another way to introduce these 
polynomials is to require that every qkix) has the Fourier transform e~* /"^Pkiit), so that Um{t) 
appears as the Fourier transform of (pm{x). 

The following statement is standard: in the interval \t\ < n^/^ 



\Vnit) - Umit)\ < 



„(m-2)/2 

where £„ do not depend on t and satisfy £n ^ 0, as n ^ 00 (with certain powers m' and m") 
Similar bounds also hold for the derivatives of orders p = I, . . . ,m, namely 



dP , , dP 
^^n(i)-^Mi) 



< 



(m-2)/2 



n 



This bound is proved in Petrov [Pel], cf. also [Pe2], pp. 209-211 (for m > 3 and \t\ < n^/^). 
We refine this result with general values of s > 2 by proving that 



dP , , dP 

-^,^n{t)-—U^{t) 



£r 



- n(^-2)/2 



(1.^ 



However, the error term in this approximation is still not sufficiently small for our applications, 
and we have to look for other related representations of Vn- In analogue with Um, introduce 

m-2 



k=l ^ 



Theorem 1.3. LetE\XY < +00 (s > 2) . For all p = 0,1 ... , m, and all \t\ < cn}/^, 



— iVn(t) — Um(t)) = n — — 

dtp ^ ^ ' ^ ^' dtp 



with 



K|<(l + |tr')e-*'/2 



C 



+ 



£r, 



iS-2 



(1.9) 



(1.10) 



^(m-l)/2 jj^c 

Here C , c and £„ are positive constants, depending on s and the distribution of X, such that 
— >■ 0, as n —>■ 00 . 
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Thus, the closeness of Cm to v near zero determines the rate of approximation of v^s by 
the functions u^'s (which have a different formal nature). This representation will be of use in 
the proof of Theorems 1.1-1.2, since, as we will see, the Liouville integrals may be applied to 
give a point wise bound on the (inverse) Fourier transforms within the class of the functions 
of the form V{^) e~*^/^, such as in (1.9). 

Note that for s > 3 the expression in the last brackets of (1.10) is dominated by Ctt.^^"*"^)/^, 
while in the range 2 < s < 3 the second summand e„n^'^*'^^) dominates the first one. In any 
case, with respect to the growing parameter n, the bound (1.9) is sharper than the one given 
in (1.8). This observation explains the improvement of (1.4) compared with the relation (1.3). 

We also remark that Theorem 1.3 holds for a more general class of functions, including 
Fouricr-Sticltjcs transforms v{t) of finite (signed) measures with finite s-th moment, such that 
v{0) = 1, v'{0) = 0, v"{0) = — 1. For example, the approximating functions Um and are 
not positive definite, but belong to this class. 

The exposition of this paper is based on chapters of auxiliary results, which are organized 
in accordance with the following table. 

Contents 

1. Introduction. 

2. Differentiability with improved remainder terms. 

3. Differentiability of Fourier- Stieltjes transforms. 

4. Cumulants. The functions tpzit) = + ^ logv{tz). 

5. The case of moments of order 2 < s < 3. 

6 Definition of the expansion polynomials Pk- 

7. Associated projection operators. 

8. Bounds of Pk and their derivatives. 

9. Edgeworth-type expansion for the functions v(tz)^^^ . 

10. Proof of Theorem 1.3. 

11. Liouville fractional integrals and derivatives. 

12. Fourier transforms and fractional derivatives. 

13. Binomial decomposition of convolutions. 

14. Proof of Theorems 1.1 and 1.2. 

2. Differentiability with Improved Remainder Terms 

For our purposes, we use the following terminology. 

Definition 2.1. Let a complex- valued function y = y{t) be defined in some interval 
a < t < b, and let s > 0. We say that y is s-times differentiable, if it has continuous 
derivatives up to order m = [s] in (a, b), and for any to £ (a, b), as i — >■ to, 

yim)^t) = y^^\to) + 0{\t - torn- (2-1) 

The case s = corresponds to continuity, while the case of a positive integer s = m - to 
the property of just having continuous derivatives up to order m. 
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The following obvious characterization will be an important tool in the derivation of the 
Edgeworth-type expansions for charactersitic functions. It is obtained from (2.1) by the re- 
peated integration over the variable t near to. 

Proposition 2.2. Let y have continuous derivatives of order up to m = [s] in (a, 6). The 
function y is s-times differentiable on {a,b), if and only if, for any point to G (a, 6) and all 
p = 0, . . . ,m, as t to, 

fe=0 

One can also provide quantitative estimates on the remainder term in (2.2), if we start with 

\i"'Ht)-y^"'\to)\<\t-tor"'ei\t-to\), 
where £ = e{w) is a non-decreasing function in u; > 0, such that e{w) — > 0, as — >■ 0. Then 



dP dP ^ yW(to) , 



<\t-tor^ei\t-to\), 



k=0 

for any p = 0, . . . ,m. 

By the chain rule given below as Lemma 2.4, we have the following: 

Proposition 2.3. If y is s-times differentiable on {a,b), s > 0, and z = z{y) is analytic 
in some domain, containing all values y{t), then the superposition z{y{t)) is also s-times 
differentiable on {a,b). 

Lemma 2.4. Under the conditions of Proposition 2.3, z{y{t)) has derivatives up to order 
m = [s] on {a,b), given by 

. . ^^ ,^d^^+-+^p z{y) ^ I f I d'y{t)\''^ , , 

- z{y{t)) =p\^ ^^^^^^ n k;\v\^) ' 

for all p = I, . . . ,m, where the summation is performed over all non-negative integer solutions 
(fci , . . . , kp) to the equation ki + 2k2 + • • • + pkp = p. 

Proof of Proposition 2.3. By the definition, for all to G (a, b) and r = 1, . . . , m, 

d'yit) 



df 



cr + oilt-tol''""), as t^to, 



where = y^'^\to)- Raising these equalities to the fc^-powers and then multiplying them, we 
get a similar representation 



^ 1 fl d^y(t)\'''^ 

n^U^j =^^^oi\t-to\^ 

r=l ^ ' 
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with some constants c^, depending on the m-tuples k = {ki, . . . , km)- In addition, putting 
j = ki + ■ ■ ■ + km, we have 

z^^\y{t)) = z^^\y{to)) + 0{\y{t)-yito)\) 

= z^^\y{to))+o{\t-tor"'). 

Inserting these relations in (2.3) with p = m, we obtain that [^(y(t))]("^) = c + o{\t - iol"""") 
with some constant c. 

In addition, the right-hand side of (2.3) represents a continuous function in t, so necessarily 
c = [z{y)]^™-'' (to). This means that (2.1) is fulfilled, and Proposition 2.3 is proved. 

3. Differentiability of Fourier- St ieltjes Transforms 

A large variety of examples of s-times differentiable functions appear as Fourier-Stieltjes 
transforms of finite measures on the real line with finite absolute s-th moment. 

Proposition 3.1. Let X be a random variable with characteristic function v{t) = Ee**"*^. 
//E < +00, s > 0, then v is s-times differentiable on the real line. Moreover, its m = [s] 
derivatives are representable, as t 0, by 

v^\t) = ^ E {iXy+'^ - + o(|ir-^), p = 0, . . . , m. (3.1) 
fe=o 

One can state a similar proposition for more general Fourier-Stieltjes transforms 

/ + 00 
e'^'^dFix), 
-oo 

where F is a function of bounded variation on the real line, such that / d |F|(a;) < +oo 

(where \F\ denotes the variation of F, viewed as a positive finite measure). On the other 
hand, such a more general statement may be obtained from Proposition 3.1, as well. Indeed, 
one can always represent F as a linear combination ciFi — C2-F2 of two orthogonal probability 
distributions (with ci,C2 > 0). Then |F| = ciFi + C2F2, so J \x\^ dFi{x) < +00. Applying 
Proposition 3.1 to Fi, we obtain a similar statement for F. 

Proof of Proposition 3.1. By the moment assumption, the characteristic function v has 
m continuous derivatives, described by 

/+00 
e'^"" {ix)PdF{x), p = 0,l,...,m, 
-00 

where F is the distribution function of X. In particular, 

/+00 
e^'^ixr dF{x). (3.2) 
-00 

Hence, the relation (3.1) would follow immediately from the formula (2.2) of Proposition 2.2, 
once it has been established that v is s-times differentiable. Namely, we need to see that, for 
any to G H, as t ^ to, 
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The formula (3.2) is telling us that v^'^^ represents the Fourier-Stieltjes transform of 
the finite measure Fm{dx) = {ix)^ dF{x) with J \xY'''^ d\Fm\{x) < +00. Representing 
Fm = ciG + C2H with distribution functions G and H, we can also represent f '-"^^ as a linear 
combination of the two Fourier-Stieltjes transforms with J \xY~'^ dG{x) < +00 and similarly 
for H. Hence, in order to show that v is s-times differentiable, it is enough to consider in (3.3) 
the case m = and < s < 1, only. 

Thus, Proposition 3.1 has been reduced to the case E|X|* < +00, < s < 1, when one 
needs to show that 

v{t) = v{to) + o{\t-tQ\'). 

Moreover, without loss of generality, it suffices to consider the point to = 0, only, in which 
case we need to show the relation v{t) = 1 + o(|t|*)- The case s = is immediate, so let s > 
and write 

/+00 
(l-e^*^)dF(x). 
-00 

For definiteness, let i > 0. Since in general |1 — e*^| < min{2, |a;|}, a; G R, we have 
\\-v{t)\ < 2 / dF{x)+t I \x\dF{x) 

J\x\>2/t J\x\<2/t 

= 2 [ dG{x)+t f xdG{x), 

Jx>2/t Jx<2/t 



where G is the distribution of |X|. By the assumption, 

r+oo 

ip{x) = B\X\'l^\x\>x}= y'dG{y)^0, as x ^ +00. 

J X 



We have -^{x) > x' JJ°° dG{y), so J^^ dG{y) = o{x-'), that is, J^^^^^ dG{x) = o{t'), as 
t ^ 0. Finally, by integration by parts. 



Hence, 



t [ xdG{x) <f f 

Jx<2/t Jo 



dy. 



ix<2/t Jo y 

But the last integral tends to zero, as t ^ 0, by the Lebesgue dominated convergence theorem. 
Proposition 3.1 is proved. 



4. Cumulants. The functions tpzit) = + p ^og v{tz) 

If a complex valued function v on the real line has m continuous derivatives and v{0) 7^ 0, 
then v{t) 7^ in some interval (— c, c). Moreover, the principal value of the logarithm logt;(t) is 
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well-defined in that interval and represents a function, which has also m continuous derivatives. 
The corresponding derivatives at the origin (with a proper normalization), 



k = 0,1, . . . ,m, 



will be called the generalized cumulants or just cumulants, associated to v. 

This terminology is standard, when v represents the characteristic function of a random 
variable (with m finite absolute moments). However, we shall have more general classes of 
functions. 

Applying Propositions 2.2-2.3, we arrive at: 

Proposition 4.1. Let v be s-times differentiable on the real line, s > 0, not vanishing in 
some interval (— c, c). Then, logv is s-times differentiable in (— c, c). In particular, as t 0, 

^logv{t) = ^Y.^-^{ttf + o{\trn, p = Q,l,...,m, m=[s]. (4.1) 

k=G 

Note that if v has m + 1 continuous derivatives, then, by the usual Taylor's theorem, the 
remainder term in (4.1) can be sharpened, and we have 

^ \ogv{t) = ^ E ^ at? + 0(|t|(-+i)-^'), p = 0, 1, . . . , m. (4.2) 

If v{t) = Ee**^ is the characteristic function of a random variable X, the assumptions of 
Proposition 4.1 are fulfilled, as long as E|X|* < -|-oo (Propositions 3.1). Then 7^ are usual 
cumulants with 70 = 0, so (4.1) should be written 

^ = ^ £ It + '^^l*!^"')' = 0' 1' • • • ' 

k=l 

In the particular case p = 0, the above equality becomes 

m 

logv{t) = Y,y_i}tf + om. (4.3) 
k=l 

However, the case p > 1 in (4.1) cannot be deduced directly from (4.3). 

Let us recall how to relate the cumulants to the moments = EX^. Applying Lemma 2.4 
with z{y) = logy, y = v{t), and at the point t = 0, one obtains a well-known identity 

7, = p! E (-1)^=^-^-+'=-^ (^1 + • • • + - n ^ (fr) 

r=l ^ ^ ' 

where the summation extends over all non-negative integer solutions {k\, . . . , kp) to the equa- 
tion ki + 2^2 + • • • + pkp = p. Note that 7^ depends on the first p moments of X, only. For 
example, 

71 = "1) 73 = as - 3aia2 + Oi, 

2 2 2 4 

72 = a2 — oil, 74 = 04 — 3a3 — 4aia3 -|- 12aiQ;2 ~ Goii- 
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Under standard moment assumptions, such as EX = 0, EX^ = 1, we have 71 = 0, 72 = 1, 
73 = CKs, 74 = ^4 — 3. For any normal random variable, 7fc = 0, for all A; > 3. 

Now, returning to the general setting as in Proposition 4.1, assume that s > 2, and v{0) = 1, 
^;'(0) = 0, v"{0) = -1. Then 70 = 71 = 0, 72 = 1, so logv{t) = + od^H- Therefore, it is 
natural to center and normalize this function by introducing the family of the functions 

i^zit) = ^t'^ + log v{tz), \tz\<c, 

where 2: 7^ is a given parameter. Clearly, ipz is s-times differentiable in this t-interval, and 

MO) = i''M = <(0) = 0, ^i^')(0) = 7p i'zP-^ {p = 3,...,m). 

Moreover, reformulating Proposition 4.1 in terms of the functions t ^ \ogv{tz) with fixed 
2; 7^ 0, we get: 

Corollary 4.2. Let v{t) he s-times differentiable on the real line, s > 2, not vanishing for 
\t\ < c, and such that t;(0) = 1, v'(fd) = 0, f"(0) = —1. Given z ^ Q, in the interval \tz\ < c, 
for all p = 0,1, ... ,m, m = [s], 

t Mt) = ^ E ^ + \tr^ i^rMtz), (4.4) 

fc=3 

where s = £{t) is defined and continuous in \t\ < c and satisfies e{t) — )■ 0, as t — )• 0. 

Also, as remarked after Proposition 4.1, cf. (4.2), when v has m + 1 continuous derivatives, 
a representation with sharper remainder term holds, 

^ ^^(t) = ^ y ^ ^it)k ,k-2 ^ |^|(m+l)-p |^|m-l 

dtp ^ ' dtp ^ k\ ^ ' ^11 II 

fe=3 

where A = A{t, z) is a bounded function in the region \tz\ < c. 

In the case of characteristic functions. Corollary 4.2 admits a slight refinement. 

Corollary 4.3. Let X he a, random variable with the characteristic function v, and with 
EX = 0, EX^ = 1, E\X\' < +00, s>2. Then, given z 0, the relation (4.4) holds in the 
interval \tz\ < \/2. 

Indeed, by Taylor's theorem, |1 — v{t)\ < \ t^, for all t G R. Hence, one may choose in 
Corollary 4.2 any value < c < \/2. 

5. The Case of Moments of Order 2 < s < 3 

In case 2 < s < 3, Corollariy 4.2 is simplified, since then there are no terms in the sum 
(4.4). In particular, when s = 2, we have 

^,{t) = t'eQ{tz), '^',{t) = \t\ei(tz), i^lit) = e2{tz) (5.1) 



Edgeworth-type expansions for charetcteristic functions 



11 



with some functions ei{z) — )• 0, as 2; — )• 0. This leads to the following observation which is 
classical in case of characteristic functions. 



Proposition 5.1. Assume thatv{t) has two continuous derivatives withv{0) = 1, v'{0) = 
and v"{0) = —1. There is a function — >■ +00, for z — >■ (0 < < 1), such that uniformly 
in the intervals \t\ < 



dtp 



v{tz[ 



dtp 



<e-''/^e{z), p = 0,1,2, 



(5.2) 



where e{z) — )■ 0, as 2; — )• 0. 



Proof. For completeness we include a well-known argument. Let \t\ < c be an interval 
where the function v{t) is not vanishing. Choose the function T = as to satisfy Tz\z\ < c, 
whenever < < 1, and Tz\z\ — >■ 0, for z — >■ 0. These conditions will be assumed from now 
on. Moreover, for any continuous function V{t), one can choose — >■ +00, for z ^ 0, such 
that 

sup ^0, as 2; ^ 0, 

\t\<T. 

and similarly for the first two derivatives of ipz. 

For the proof, it is enough to see that, whenever e{z) and W{t) > is continuous 
and increasing in t > 0, one can choose Tj — +00, such that W{Tz) sup|j|<2^^ ~^ 0, as 

z — 7- 0. Here, we may assume in the following that e{z) > is even and also increasing in 
2; > 0. Then, the latter statement may be simplified to WiT^) e{Tzz) — )• 0, which is obviously 
true with a sufficiently slowly growing T^. 

In particular, in view of (5.1), with some +00, as z ^ {0 < \z\ < 1), we have 

e{z) = sup (IV'.(t)l + WM + WM) ^0, as z ^ 0. (5.3) 

\t\<T. 

Now, write v{tzy/^'^ = g{t)e'^''^^\ where g{t) = e~*^/^. Applying (5.3), we get 

\v{tz)'/'' - git)\ < git) le'^^W - 1| < Cg{t)\Mt)\, 

for all \t\ < Tz with some constant C. Since also ipz{t) ~^ uniformly in that interval, we 
arrive at the desired conclusion in case p = 0. 

Writing {v{tzYl^ )' = g' {t)e'^'^^^^ + g{t)ip'^{t)e'^'''^*^ and using the previous step, wc get 

\{v{tz)y^'y - g'{t)\ < C\t\g{t) \Mt)\ + Cg{t) \^'M- 

Since ipzi't) ~^ ^^'^ tipzit) — >■ uniformly in that interval with an appropriate choice of T^, 
we arrive at the conclusion in case p = 1. The case p = 2 is similar. 

Now, let us turn to the range 2 < s < 3. In this case, we obtain up to polynomial factors 
m front of e-*''^ in (5.2) more information about the possible growth of T^. 
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Proposition 5.2. Let v{t) be s-tirnes diff'erentiable, 2 < s < 3, not vanishing for \t\ < c 
(c > 0), and such that v{0) = 1, v'{0) = 0, v"{0) = -1. Given < |z| < 1, 



— vitz)'-''' - -;—e ' 
dtp ^ ' dtp 



< {\t\'-^ + \t\'+^)e'' /^s{z), p = 0,l,2, 
uniformly for \t\ < c \z\''^^~^^/^ with some function e{z) — >■ 0, as z — >■ 0. 



Proof. By Corollary 4.2, in the intervals |t| < = c \z\ ^^/^ with < < 1 

\Mt)\ < \t\' \zr^ e{z), 

< \tr'\zr'e{z), 
mt)\ < \tr'\zr'eiz) 

with some bounded function e{z) — t- 0, as z — t- 0. Indeed, the conditions < jz] < 1 and 
\t\ < Tz insure that \tz\ < c|2:p/* < c and also \tz\ — )• 0, as z — )• 0, uniformly in \t\ < Tz- 
Now, by the first inequality, 

\^l^zm<W\zY-^e{z)<C {\t\<Tz) 

with some constant C. Hence, using the same notations and arguments as in the proof of 
Proposition 5.1, for all |t| < T^, 

\v{tzfl'" -g{t)\ < Cgm^m < C'g{t) ■ \tY \zY-'' e{z). 
Since also IV'UOl — with some constants C, C" we get 

\{v{tzf'^')'-g'{t)\ < C\t\g{t)\Mt)\+C9{t)\^'z{t)\ 

< C'{\t\'+' + \tr') git)eiz). 

Finally, writing 

{v{tz)'/^')" = g"{t)e^^('^ + 2<7'(t)<(t)e^^W + 2git){^': + ^p'^itf) e'^^W 
and using IV'zlOl ^ \^\^~'^ "^^ that up to some constants 

{v{tz)'/^y - g"{t)\ < Ct'g{t) \Mt)\ + C\t\g{t) + Cg{t) {W^{t)\ + I^UOH 

< C ^ ^ |^|s-2 ^ |i|2(s-l)^ \z\'-^e{z). 

All powers of \t\ vary from s — 2 to s + 2, so Proposition 5.2 is proved. 



6. Definition of the Expansion Polynomials 

The polynomials Pfc introduced in Section 1 appear not only in connection with character- 
sitic functions, but in a more general setting as well. 

Namely, let v{t) be a complex- valued function on the real line, which is s-times differentiable 
(s > 2) and such that f (0) = 1, u'(0) = 0, u"(0) = — 1. Then v has cumulants 7^, = 1, . . . , m, 
where m = [s]. Moreover, 71 = and 72 = 1. 
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Assume v is not vanishing in the interval \t\ < c, and let us return to the functions 



where z 7^ is viewed as a (small) parameter and \tz\ < c. Recall that, by Corollary 4.2, 

e*'/2 v(tzy/^' = eMt) = exp i J] z'' + \t\^zr'' e{tz) \ , (6.1) 

I- k=i (^ + 2)- J 



where E{t) is defined and continuous in \t\ < c and satisfies e{t) ^ 0, as i ^ 0. Moreover, if 
V has m + 1 derivatives, the remainder term here may be replaced with 
where A = A{t, z) is bounded in the region \tz\ < c. 

The sum in (6.1) is vanishing in case 2 < s < 3. To study this representation in the case 
s > 3, introduce the polynomials 

m-2 

By a formal Taylor's representation with respect to the (complex) variable z, we have 

00 

e^'^W = l + ^afe(it)/ (6.2) 
k=l 

with some coefficients ak{it). To justify this step and precisely determine the coefficients, write 



00 



E 



I ■ 



p=0 P 



which can further be expanded as 



^ Pi! . . . Pm-2! • V • ; 

p=0 Pl+-+Pm-2=P ^ 

The whole double sum is absolutely summable for all complex numbers t and z. Indeed, let 
C = Y1T=3 • fixed integer p > 0, the finite sum of the absolute values of the terms 

in (6.3) is bounded by 



M y ll^l^l^^+si^lfc] <icM max \t\^+^\zA . 
p!V^(fe + 2)!'' / ~ p! \i<k<m-2^ ^ ''y 



Assume without loss of generality (in order to get some quantitative bounds) that \t^z\ < 1 
and l^l < 1. Then, 

|^|fc+2 i^jfc < |^|-(fc+2)/3 |^|fe ^ |^|(2fc-2)/3 < 

Hence, the above sum is bounded by C^/p\ Furthermore, note that |W2(t)| < C. 

Thus, the total sum of the absolute values is bounded by e'^ , and one may freely choose 
the order of summation. Collecting the coefficients in (6.3) in front of z'^, we arrive at (6.2) 
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with 



1 / \ pi / \Pm-2 



Pi+2p2H h(m-2)pm-2=fc 

where the summation is extended over all non-negative integer solutions . . . ,Pm-2) to the 
equation pi + 2p2 + • • • + (m — 2)p^_2 = 'i:- Note that 

3pi + 4p2 H h mPm-2 = k + 2(pi + P2 H h Pm-2)- 

Hence, replacing with t, 

pi+2p2H \-{m-2)pm-2=k 

In addition, ii k < m — 2, the condition pi + 2p2 + • • • + (m — 2)]3^_2 = k implies that Pk+i = 
■ ■ ■ = Pm-2 = 0. Therefore, in this case depends on the first k cumulants 73, ... , 7^+2^ only. 
More precisely. 

For example, if m = 3, we have ai{it) = ^ {it)^- In case m = 4, ai(i) = ^ and 

^ ' 72 24 

In general, subject to p\ + 2^2 + • • • + kp]^ = k, the expression k + 2{pi + • • • + pk) does 
not exceed 3A; and reaches this value (when pi = k, p2 = ■■■= Pk = 0), so ak represents a 
polynomial in t of degree exactly Sk. 

Definition 6.1. Given an integer m > 3 and complex numbers 73, . . . ,7m, one defines 
(1 < A; < m — 2) as the polynomial introduced above, namely, 

" ^ pi!...pfc!U!y "\{k + 2)\) ^ 

Pi+2p2+-+fcpfe=fc'^^ K/c \ / \V I / / 

With this definition, the representation (6.2) may also be written as 

f m -. m—2 00 

exp ^ ^ {itfz^-^ = 1 + E + E «^(^*)^'' (6.5) 

fe=3 ■ ^ k=l k=m-l 

where a^'s are described in (6.4). 



7. Associated Projection Operators 

Let us note first that every polynomial in (6.4) contains terms involving powers of t not 
smaller than k + 2 (since necessarily pi + ■ ■ ■ + Pm-2 > !)• This observation may be used to 
obtain an initial trivial bound for the last sum in (6.5) in case of small values of t. 
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Lemma 7.1. Given complex numbers 73, . . . ,7m) for all complex z and t, \t\ < 1, 

00 

J2 Wk(it)z''\ < (e^^^^ - 1) \tr^\ (7.1) 

k=m—l 

Where C{z) = ^]:^,\^\z\>'-\ 

Indeed, using (6.4) and \t\ < 1, we see that each term in the sum of (7.1) is bounded by 
up to the factor 



E 



Pi /U, I U|m-2\ Pm-2 



ITsII^IV /^|7m|k 



_L 9^ . Pi!... Pm-2! V 3! y \ ml 

pi+2p2H |-(m-2)p,„_2=K 

Summation of these expressions over all > 1 results in e'-'^^^ — 1. 

The bound of Lemma 7.1 is needed in order to express the "cumulants" jk directly in terms 
of the associated polynomials Pj. 



Lemma 7.2. Given complex numbers 73, . . . , 7m, for all complex z and k = 3, . . . ,m, 

(7.2) 



k\ i^ dt^ 



.7=1 ^ 



t=0 



To see this, note that, by (6.5) and (7.1), we obtain that, as t ^ 0, 

m /■ m—2 \ 

^ ^ ^u)'^,k-2 ^ Pjiit)z^j + o{\tr+'). 

k=3 j—^ 

Since the right-hand side (without the remainder term) represents an analytic function near 
zero, comparison of coefficients powers of t immediately leads to (7.2). 

The identity (6.5) suggests to consider special operators, defined on the space Vg of all 
s-times differentiable functions u : R ^ C, s > 2, such that v{0) = 1, ^'(O) = 0, v"{0) = —1. 

Definition 7.3. Given v £ Vs, s > 2, and an integer 2 < m < s, we put 

m-2 



( III — ^ \ 
l+J2Pk{it)], teR, 
k=i ' 



where are the polynomials from Definition 6.1 based on the cumulants 

Tfe = ^fe^ log^(*)lt=o' A; = 3,..., m-2. 

If m = 2, there are no cumulants and polynomials in the definition, so {T2v){t) = e~*^/^, 
for any v &Vs. If m = 3, 



16 S. G. Bobkov, G. P. Chistyakov and F. Gotze 

for any v £ Vs, s > 3 (where 73 may be an arbitrary complex number). If m = 4, then for any 
veVs,s>A, 

(T.vm = e-^y^ (1 + 1 i^tf + I i^tr + | i^tf^ . 

Clearly, every TmV is an entire function and hence belongs to all V^, s > m. This defines 
an operator Tm : Vg ^ Vg which turns out a projection operator 

Proposition 7.4. We have T„iT,nV = T„,,v for any v G Vs, 2 < m < s. Moreover, T^v 
and V have identical derivatives at the origin up to order m. 

Proof. The statement is equivalent to the property that = T^u and v have equal 
cumulants. Let 7^ and 7^ denote the cumulants of and respectively (3 < A; < m). By 
Definition 7.3, 

But the right-hand side equals according to Lemma 7.2 applied with z = 1. 
Thus Proposition 7.4 is proved. 

Note that T^v does not need to be a characteristic function, even if u is a characteristic 
function of some random variable. However, it always represents the Fourier-Stieltjes transform 
of a finite signed measure. 

In the following we approximate v by its projections T^f. Combining Proposition 7.4 with 
Proposition 2.2, we get: 

Corollary 7.5. Given v eVg, 2 < m < s, as t ^ 0, 

— [v{t) - Tmv{t)) = o(|ir-^), p = 0, 1, . . . , m. 

Finally, let us formulate an asymptotic property of the projection operators for growing 
parameter m (although this will not be needed in the sequel). 

Proposition 7.6. Assume that v{t) admits an analytic extension to the disc \t\ < p, where 
it has no zeros, and v{0) = 1, v'{0) = 0, v"{0) = —1. Then T^v{t) — >■ v{t), as m 00, i.e., 

vit) = e-''/^(l + f2Pkiit)\, \t\<p. 
^ k=l ^ 

Moreover, the series is convergent absolutely. 

If v{t) = Ee**"'^ is the characteristic function of a random variable X, the assumptions of 
Proposition 7.6 are fulfilled, provided that EX = 0, EX^ = 1, Ee''l"^l < +00 (that is, an 
exponential moment of order p is finite) and v{t) does not vanish in the disc \t\ < p. 
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Proof. By the assumption, logf(t) is analytic in the disc \t\ < p, so it is representable as 
the sum of the absolutely convergent power series 

oo 

logv{t) = Y,^iit)\ \t\<p. (7.3) 

k=3 

Hence, starting with (6.5) with z = 1 and letting there m — >■ oo, it is sufficient to show that 
Ylk^m-i l^^kiit)] — >■ (note that a^'s also depend on m). 
Rewrite the representation (6.4) as 



pi+2p2H 'r{m-2)pm-2=k ' 

which implies that 

Pi+2p2H \-{m-2)pm-2=k 

Here the right-hand side may be bounded by the quantity 



Pl+2p2+3p3+---=fc r=l 



Pr\ V (r + 2)! 



which does not depend on m. After summation over all > 1 (thus removing any constraint 
on pr), we get h{t) = e^d*!) - 1, where C{a) = ^^^3 ^ aK But C{\t\) < +00 for all 

\t\ < p in view of the absolute convergence of the series (7.3). Hence, in this case 

00 00 

\ak{it)\ < ^ bk{t) ^ 0, as m — >■ 00. 

fc=m— 1 fc=m— 1 

With similar arguments, we also obtain that Ylh=i \-^k{it)\ < +00 for \t\ < p, in view of 
Definition 6.1. Thus Proposition 7.6 is proved. 



8. Bounds of Pk and their Derivatives 

We will need a bound similar to the one in Lemma 7.1 which extends to large values of t 
and involving derivatives of the polynomials and Pk- 

To this aim, we start with arbitrary complex numbers 73, ... ,7^, m > 3 (which may be 
interpreted as cumulants of a given function v) and return to the representation (6.2), 

00 

u;^(i) = e^-W = l + ^afc(^^)^^ t,zeC, (8.1) 

k=i 

where 

m-2 

^^W = E(fc^(^^)'^'^'' (8-2) 



18 S. G. Bobkov, G. P. Chistyakov and F. Gotze 

and where the polynomials are described in (6.4). By the very definition, ak = Pk, as long 
as A; < m — 2. 

As we have already noticed, the sum in (8.1) is absolutely convergent and therefore repre- 
sents an entire function with respect to z, for any fixed t. It is also clear that the series may 
be term wise differentiated, so that 



w(P\t) = Y,i''a^k\it)^'^ p>l, (8.3) 
k=l 



which is absolutely convergent, as well. 

In order to bound and its derivatives, we use the quantity 

m 
k=3 

One natural approach (which is however different than the one in [Pe2]) is based on the 
application of Cauchy's integral formula 

with a suitably chosen parameter p > 0. In view of (8.3), there is a more general identity, 
involving the derivatives. 

Lemma 8.1. Let \tz\ < 2 and \t^z\ < 2, and let < p < m be an integer. Then 

\W^\t)\ < 2"*-2C|ir^'min{l, \t\'^}. (8.5) 

Indeed, by definition (8.2), 

w^p\t)= y -^ii-i^-fz''-2. 

q=m.ax.(p,S) 

But \ti-Pzi-'^\ = \tz\i-^\t^z\\t\-P < 2"'-^\t\-P, whenever 3 < q < m. Hence, \W^\t)\ < 
2m-2(-j On the other hand, just using \z\ < 2/\t\, we get 

29-2 

U9-P^9-2| < Uiq-P^_^^ 2«-2|t|2-P < 2"^-^\t\'^-P . 

This gives an improvement over the previous estimate in case \t\ < 1 and proves (8.5). 
Lemma 8.2. For all integers k > 1, < p < m, and all complex t, 

\afm < Cm,,\t\-p min{i,m (^ ^--m^m y 

with constants Cm,p = (4"'(1 + C))p e^"^ . 
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Proof. Given t / 0, we choose in (8.4) the radius 

2 

^ ~ max{|t|, |t|3}' 

Hence, on the circle \z\ = p, both \tz\ < 2 and \t^z\ < 2 arc fulfilled, thus the inequality (8.5) 
may be applied. In particular, |W2(i)| < 2'"~^C, and from (8.4) with p = we get the desired 
estimate 

\ak{it)\ < ^e^"*"'^. 
Next, by the formula (2.3) of Lemma 2.4, for all p > 1, 

where the summation is taken over all non- negative integer solutions {ki, . . . , kp) to the equa- 
tion fci + 2k2 + • • • + pkp = p. Hence, using (8.5), given that \z\ = p, we arrive at 

< e^-'^ir-P»En^( '"'"^7^^''"^' )^- 

r=l ^ ^ ' 

Since necessarily 1 < A;i + • • • + pkp < p, the product may be bounded by the product of 
min{l, (in the first power) and 2"*~^C (replaced by 2"*~^(1 + C)), and raised to power p. 
This leads to 

\w^Ht)\ < e^""^^ {2'^-'^{l + C)Y\t\-Pm.m{l,\t\'^}Bp, (8.6) 

where Bp = P^-J2Tlr=i ^(h)'^'^- This constant can also be described by virtue of the same 
formula (2.3), applied with z = and y{s) = e^, in which case it reads 

r=l ^ 

One should apply this formula at ,s = 0, thus we consider the functions 6p(s) = (e^°)*-^\ 
s > 0, and their values hy = bp{0) = Bp/e. The recursive identity 6p+i(s) = (e^^e^^)*^^^ = 
J2r=o ^r{s) implies that the sequence r — )• 6^ is non-decreasing and 6p+i < 2Pbp. There- 
fore, 

bp < 2P-\^i < 2P-^2P-\^2 <■ ■< 2P-^2P-^ ...2°bo = 2P^-'^y^ e. 
Hence, Bp < 2p(?'-i)/2, and together with (8.6) this gives the estimate 

\wiP\t)\ < e2'""'^(2(— 2)+(f-i)/2(i + c))^'|i|-P^in{i,|i|2|. 
It remains to apply (8.4) and simplify the constant. Thus Lemma 8.2 is proved. 
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Now, fix an integer p = 0,1 ... ,771, and assume that < f = max{|t| \t\^} ' ^^^^ l^'^l — ^ 
and \t^z\ < 1. By Lemma 8.2, 

oo oo I 1^ 

\a'^^\it) z'^l < Cm,p\trrain{l,\t\'} J] ^ 

k=m—l k=m—l 

< 2C^,p|t|-fmin{l,|t|2} 



pm-l 



2; 



< C^,p|t|-^'min{l,|t|2} (max{|t|,|t|=^}) 
To simplify the dependence in t, note that in case \t\ < 1, 

\t\-Pmm{l,\t\^} (max{|t|,|f|3})'""^ = 

while the left expression is equal to |i|3(™-i)-P in case \t\ > 1. 

Also note that the condition \tz\ < 1 is fulfilled automatically, as long as \t^z\ < 1 and 
\z\ < 1. Therefore, recalling also that Pk = Uk ior k < m — 2, we obtain: 

Proposition 8.3. // < |z| < 1 and \t^z\ < 1, then for all p = 0,1 ... ,m, 

^ k=l ' 

where \A\ < Cm,p = (4"*(1 + C))p e"^""^ . 

9. Edgeworth-type expansion for the functions v{tzY/^^ 

Assume that v{t) is s-times differentiable, s > 2, and not vanishing for \t\ < c (c > 0), and 
such that v{0) = 1, v'{0) = 0, v"{0) = -1. For 

vS) = v{tzf''" 

define the approximating functions 

m-2 

UrnXt) = Um(t,z) = e " I 1 + i^knt)z- 1 , m= \s\ 



(iri—^ \ 
1 + J]p,(zt)zM, 
k=i ' 



where the polynomials arc based on the cumulants 73, ... , 7^ of v. Put m'(s) = s — p, 

m"{s) = 3(m — 2) + max{s + p,{s — l)p}. 

In particular, m'{s) = s and m"{s) = s + 3(m — 2) in case p = 0. Note that m"{s) < 2m^ in 
all admissible cases. 

In this section, the relation (1.8) is established in the following more general form. 

Proposition 9.1. Let s > 3. Given z real, < < 1, in the interval \t^z\ < c^, for all 
p = 0,l,...,m, 

\v^\t) - u^P\t)\ < {\tr' + \tr") e-*V2 \zr^e{z), (9.1) 
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where e{z) — )■ 0, as z — )■ 0. Moreover, if s > 2 and v{t) has (m + 1) continuous derivatives, 
then with some constant A, and with m' , m" , corresponding to s = m + 1, 

\v^\t) - u^\t)\ < A {\tr' + \tr") e-*'/2 |^|m-i_ (9 2) 

We will refer to (9.1) and (9.2) as the scenarios 1 and 2, respectively. Note that in the 
second case, although v has cumulants up to order m + 1, we require that 7m+i does not 
participate in the definition of the polynomials Pfe. In particular, the value m = 2 is covered 
in (9.2), and then Um{t) = e~* (that is. Pi is not present). 

Proof. Without loss of generality, assume c = 1. Write Vz{t) = e~*^/^ w^z(i) e'^^(*\ where 

m 

w,it) = e^^('\ W^,(i) = ^ (^^)'^'=-^ 

A;=3 

Mt) = ^ + ^ log vitz) = log {e^'/^v,{t)) , h,it) = i),{t) - W,{t). 
By the definition of and Pfc, 

m— 2 oo 

w;^(t) = 1+ ^ Pfe(ii)z'= + P^(i), Rz{t)= ^ ak{it)z^. 

h=\ k=m—l 

Therefore, 

vAt) = u„Xt)e^^^'^ +RS)9{t)e^''''\ git) = e-*'/'. 
Given p = 0,1 . . . ,m, wc differentiate this representation according to the Leibnitz rule: 

4'\t) - u^P\t) = Ii+I, + I, = u(P\t) (e'^^W - 1) 



+ ^C^ti(^-*^)(i)(e'*^W)(*^) 

k=l 

+ j^C^v {Rz{t)9{t)f^ (e'^^(*))M, (9.3) 

fe=0 

where Cp = i^\(jLky are the combinatorial coefficients. Note that when p = 0, the second term 
I2 is vanishing. 

Estimation of Ii. 

In Corollary 4.2 it is shown that, if \z\ < 1 and \t^z\ < 1, the functions hz and their 
derivatives are uniformly bounded and admit the bounds 

\h^f\t)\ < \trP\zr%iz), p = 0, 1, . . . ,m, (9.4) 

where each ep{z) is defined in \z\ < 1 and satisfies ep{z) — )• 0, as z ^ 0. Moreover, if v has 
m + 1 continuous derivatives, then we have a sharper estimate 

|^^)(i)| < Ap l^r-i (9.5) 

with some constants Ap. In particular, when p = 0, these bounds correspondingly give 

\eh^{t) _ i| < \t\'\z\'-^e{z), le'^^W - 1| < ^0 I^P^^ \zr~\ (9.6) 
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with some e{z) — )■ 0, as z — )• 0, and a constant Aq. 

On the other hand, since every has degree 3k < 3(m — 2), and < 1, for all p = 
0, 1, . . . , m, m > 3, 

m-2 



^ (l + ^ PU^t) z'^) < C{1 + |t|3(™-2)-P) 
^ k=l ^ 



with some constant C, depending on m, p and the cumulants 7fe's. Since also 



"1 g-t2/2 

dtp 



<Cp(l + |^^e-*'/^ (9.7) 



we get, by the Leibnitz rule, 

ki^HOI <C{1 + |i|3('"-2)+P) e-*'/2^ (9.8) 

where we allow the constants depend on m, p and the cumulants 73, . . . ,7m- For m = 2, 
Um{t) = e~* so (9.8) holds in this case as well {p = 0, 1, 2). Combining this with (9.6), we 
correspondingly arrive at 

< (|tr + e-*'/2 |z|'*-2e(z), (9.9) 

with some constant A and e{z) — )• 0, as 2; — )• 0. As a result, we obtain the desired bounds on 
the first term Ii in (9.3) for both scenarios. 

Estimation of l2- 

To treat the second term, more precisely - the products Um~^\t) {e^''^*))'^^\ assume that 
p > 1. By formula (2.3), for any k = 1, . . . ,p, 

1 (h^i\t)y^ 



(g?..(t))W=e'^.Wfc!^J]±(^/ 

where the summation is performed over all non- negative integer solutions (pi, . . . ,pk) to the 
equation pi + 2p2 + • • • + kpk = k. From (9.4)-(9.5) we get for the two scenarios 

\hP{t)\P^ < \4'-''^P-^\z\^'-^^P^ Erizf-^ , \hP{t)\P^ < A^r |t|((™+i)-OPr |^|(™-l)Pr_ 

After multiplication of these inequalities over all r = 1, . . . ,/j (separately in both scenarios), 
using 1 < pi + ■ ■ ■ + Pk < k together with 

k k 
s — k < (s — r)pr < {s — l)k, (m + l) — k < ((m + 1) — r)pr < mk, 

r=l r=l 

we obtain that 

|(e^^W)W| < (Itl'^-fc + \t\^'-^')'')\z\'-h{z), (9.11) 

|(e^.W){fc)| < A(|t|(™+i)-fc + (9.12) 

with some constant A and e{z) — > 0, as 2; — >■ 0. One may combine these bounds with (9.8), 
which immediately yields 

|tx(^-*^)(t)(e'*^(*))«| < (|t|«-fe + |f|(«-l)fe+{P-fe)+3{m-2)-) g-tV2 1^1^-2 £(^), 
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Since k varies from 1 to p, the right-hand sides can be made independent of k, and we arrive 
at the desired bounds on the second term I2 in (9.3), needed for the values p > 1: 

I/2I < (|tr^ + e-*V2 \z\'-'^e{z), (9.13) 

I/2I < {\t\^"'+^^-P + |tpP+3(m-2)~j g-tV2 |^|m-l_ (-9_;^4) 

Estimation of /a. 

Now, let us turn to the third term, i.e., to the products {Rz{t) g(t))^''^ (e'*^ (*))<>-'=). By 
Proposition 8.3, for all p = 0, 1, . . . , m, 

\R'fHt)\ < C7(|t|(™+i)-P + |t|3(™-i)-P)|zr-^ 
Using (9.7) and the Leibnitz formula, the latter gives, for all k = 0, . . . ,p, 

\{R,{t)g{t)Y''^\ < C(|i|("*+^)-^ + m-'^) e-*'/2 \zr-\ (9.15) 
Case p = 0. Then necessarily k = 0, and the above inequality yields 

Ihl < C{\tr+'^ + |t|3("^-^)) e-*'/2 \zr-\ (9.16) 

It has only to be compared with (9.9)-(9.10). In the second scenario, one clearly obtains from 
(9.10) and (9.16) that 

m+l _|_ |^|5m— 5"^ 1 im— 1 



\v,{t) - Umit)\ <\h\ + Ihl < C(|tr+1 + |t|5™-5) 



This proves (9.2) in case p = 0. 

In the first scenario, in (9.16) just write = l^l''"^ £(2;) with e{z) — )• 0, as 2; — >■ 0. 

Together with (9.9) this leads to a similar estimate 

\vz{t) - um{t)\ < + l/al < C(|ir + itr+^("^-2)) e-*'/2 \zr\ 

proving (9.1) in case p = 0. Thus, Proposition 9.1 is proved in this case. 

Case 1 < p < m. li k = p, the absolute value of the product 

{R,{t)g{t)f^ (e'^^W)(P-fc) = (E,(t)5(t))(p) e'^^W 
may be estimated according to (9.15) by 

^^|^|(m+l)-p ^ |^|3(m-l)+p^ g-tV2 (9.I7) 

As in the previous step, may be replaced here with \z\^~'^ £(2;). 

If < A; < p — 1, by (9.11)-(9.12) for the two scenarios we have 

|(e?^z(t))(p-fc)| < (|t|^-(p-fc) + |t|(«-i)(f-^))|z|^-2£(2), (9.18) 

^(^ghzit)^^(j>-k)^ ^ ^ ^|^|(m+l)-(p-fe) _|_ ItpCP-'^)) l^l™'"^. (9.19) 

It remains to multiply these inequalities by (9.15). In this step we consider the two scenarios 
separately. 
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Scenario 1 (The inequality (9.1)): When multiplying (9.15) by (9.18) and looking for the 
maximal power of \t\, notice that 

{3{m - 1) + k) + {s - l){p - k) 

is maximized for k = 0, and for this value it is equal to 3(m — 1) + (s — l)p. Therefore, 

|(i?,(t)<7(t))^^'^ (e^-(t))(P-fe)| < (|t|(m+l)+{.-p) ^|^|3(m-l)+{.-l)p>)g-tV2 \z\(m-l)+(s-2)^(^^^_ 

Now, comparing with (9.9), (9.13) and (9.17), we see that the smallest power of |t| in these 
inequalities is m' = s — p. Hence, we do not loose much by writing 

\{R,it) g{t)Y''^ (e/^4t)){P-fe)| < + |i|3(™-i)+(.-i)p) ^-tV2 |^|H-i)+(.-2)^(^)^ 

which holds for all k = 0, 1, . . . ,p. To simplify, let us note that lipC"!^-!) < 1, which 
leads to 

I {R,{t) g{t)) (e'*^ W)(f-*^) I < {\tr' + |t|(^-i)P) e-*'/2 \zr'^e{z). (9.20) 
In addition, the largest power of \t\ in (9.9), (9.13), (9.17) and (9.20) is 

m" = 3{m — 2) + max{s + p,{s — l)p}. 
Hence, 73 < + e~*^/^ |2;|*~^£(2;) and similarly for Ii and I2. This proves (9.1). 

Scenario 2 (The inequality (9.2)): This case can be dealt with along the lines of scenario 1 
by letting s — > m + 1. Or, repeating the previous arguments, note that when multiplying 
(9.15) by (9.19), the expression 3(m — 1) + k + m{p — k) is maximized for A; = 0, and for this 
value it is equal to 3(m — 1) + mp. Therefore, 

\{RS)g{t)f^^ (e'^^(*))(P-'=)| < A (|t|2('»+l)-P + |i|3(m-l)+mp) ^-^ /2 |^|2(m-l)_ 

In (9.10), (9.14) and (9.17) the smallest power of \t\ is m' = {m + 1) — p. Hence, 

\{R,{t) git)^''^ (e'^^(*))(f-fe)| < A {\tr' + |t|3(m-l)+mp) g-tV2 |^|2(m-l)_ 

Again using l^l™-"^ < 1, we get 

\{R,{t) g{t))^''^ (e^^(*))(p-fc)| < (|tp' + |t|™P) e-*'/2 (9.2I) 

In addition, the largest power of |t| in (9.10), (9.14), (9.17) and (9.21) is 

m" = 3{m — 2) + max{m + p + I, mp}. 

Hence, < v4(|t|™' + |t|™") e~*^/^ and similarly for Ii and I2 in case p> I. 

This proves (9.2) and Proposition 9.1. 

One may unite Proposition 9.1 (first part) with Propositions 5.1-5.2 for the case 2 < s < 3, 
if we do not care about polynomial factors in front of 

Corollary 9.2. There is a function Tz +00, as z ^ {0 < \z\ < I), such that in the 
interval \t\ < T^, for all p = 0,1, ... ,m, m = [s], 

\v^Ht) - u'^Ht)] < s{z)\zr' e-^y^ 
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with e{z) — )■ 0. Moreover, up to some constant c > 0, one can choose = c\z\~^^^ in case 
s>3 andT^ = c \z\~^^-'^^/^ in case 2 < s < 3. 



10. Proof of Theorem 1.3 

Again, let v{t) be s-times differentiable, s > 2, such that v{0) = 1, v'{0) = 0, v"(0) = — 1. 
Note that v is not vanishing in some interval, containing the origin. 
Let us return to the family of the functions 

^ m—2 

Ur, 



(■flL — Z \ 
i + ^Pfc(ii)zM, 
fe=i ^ 



where m = [s], and the polynomials Pk are based on the cumulants 73, . . . , 7^ of w. In order 
to approximate the powers Vn{t) = one uses the values z = l/-\/n, leading to the the 

approximating functions 



(m—2 \ 
k=l ' 

On the other hand, when z = 1, we deal with the projection operators T^, i.e., with the 
functions 

(m—2 \ 
1+ ^Pfc(it)). 
fe=l ' 

Theorem 1.3 is a particular case of the following more general proposition. 
Proposition 10.1. For all p = 0,1 ... ,m, and all \t\ < cn^^^ , 



dtp ^ ^ ' ^ dtp 



+ rn (10.1) 



with 



\rn\ < (1 + 1^1'"^') e-*'/2 ( + ^1 . 

I "I - V II / \^^(m-l)/2 jiS-2 J 



(10.2) 



Here C , c and En are some positive constants, such that £„ — >■ 0, as — > 00. 



It is worthwile noting that Proposition 9.1 and thus relation (1.8) can be obtained on the 
basis of (10.1)-(10.2) as well, using the property that v{t) and e„i{t) have equal derivatives up 
to order m and both are s-times differentiable. However, we have chosen a different road of 
proof and will derive Proposition 10.1 by virtue of Proposition 9.1 (its second part). 
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In order to show how it apphes, apply the binomial formula to obtain that 



Vnit) - Um{t) = (Jnl + Crn2 + CFnS = em{—^)'^-Um{t) 



+ n 



( t \n—l 



k=2 



Thus, (T„2 is almost the term which appears on the right-hand side of (10.1), provided that 
emi^^)"'~^ is replaced with the characteristic function of the standard normal law. 

The first term Uni = ^mi:^)"' — Um{t) is of the same nature as Vn{t) — Um{t), assuming 
that Cm plays the role of v. At this point, let us recall that, by Proposition 7.4, Tjnem = em, 
and moreover, that Cm generates the same polynomials as v. Hence, Proposition 9.1, 
being applied to Cm in place of v with z = l/^/n, provides the bound on the derivatives of 
Cjni^)"' — Um{t). Since is analytic, the second assertion (9.2) of Proposition 9.1 is more 
accurate. Namely, if ejn{t) is not vanishing in the interval \t\ < c (which is true with some 
constant c > 0, depending on the cumulants, only), it gives: 



Lemma 10.2. For all p = 0,1, ... ,m and all \t\ < cn}/^ , 
where c and A are some positive constants, depending on the cumulants 73, . . . ,7m- 



(10.3) 



Moreover, using a similar argument one may estimate the derivatives of the functions 
em{t)^, which appear both in cr„2 and (7„3. Apply (9.2) with v = e^n and z = l/s/k to get 

m-2 



where A = Ak{t) is a bounded quantity in the interval \t\ < ck^^^. Putting a 
replacing the variable t with at, we obtain 



dP 



dtp 



m— 2 



+ S(l + |tp-')e-"'*'/2^ 



(10.4) 



- and 

n 



(10.5) 



where now B = Bk{t) is bounded in |t| < cn^^^ . Every Pj is a polynomial of degree at most 
3j < 3m, so all its derivatives of order up to m can be bounded by C(l + \t\Y"^ on the whole 
real line. Hence, using a < 1 and A; > 1, the same is true for the polynomial in the large 
bracket of (10.5). Using the Leibnitz rule, it then follows from (10.5) that: 



Lemma 10.3. For all p = 0,1, ... ,m, k = 0,1, . . . ,n, in the interval \t\ < cn 



1/6 



dtp "^^v^^ 



<C(l + |t|2"^ )e 



-kt^/{2n) 
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with some positive constants c and C. 
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The particular case k = n — 1 in (10.4) should be investigated in more detail by replacing 
(10.5) with a more accurate relation, namely 



dtp "^^v^^ dtp V V^i^ tt - 1)^^-'^/' 



V?^ — 1 



(assuming n > 2). Repeating the same argument concerning the growth of the polynomials 
Pj and its derivatives, and noting that, for a = we have 



dtp dtp 



C 



<hl(i + |t|P+2)e-tV2^ |t|< V^, 



n 



we arrive at the following bound (which also holds in the missing case n = 1): 
Lemma 10.4. For all p = 0,1, ... ,m, in the interval \t\ < cn}/^ 



dip 



dtp 

with some positive constants c and C. 



C 



< J_(l + |t|2™^) 

■s/n 



Finally, let us bound the derivatives of y{t) = v(t) — emit) and of its powers, which appear 
both in (T„2 and (7„3 as well. To this aim, one may appeal to Corollary 7.5, giving, as t — > 0, 



yW(t) = o(|t|^-'-), for any r = 0, . . . 



m. 



(10.6) 



In particular, y{t)^ = odt]^'^), for any k > 1. If p > 1, by the chain rule (cf. (2.3)), the p-th 
derivative of y(t)'^ represents a linear combination of the terms 

over all integer tuples (/ci, . . . , kp), such that ki + 2k2 + • • • + pkp = p and k\ + ■ ■ ■ + kp < k 
{kj > 0). By (10.6), we have b{t) = oi\t\'^), where 



S = s(^k — {ki + ■ ■ ■ + kp)) + ^^(s — r)kr = sk — p. 

r=l 

Hence, ^ {v{t) — em{t))'' = o{\t\^'^~P). Since sk — p < sm iov 1 < k < m, we obtain: 

Lemma 10.5. Let < a < ^ and c> be given. For some En 0, for all p = 0,1, ... ,m 
and k = 1, . . . ,m, we have, uniformly in the interval \t\ < cn", 



dip 



,(_L)_,„(_L))'|<,„(i + |,r)„-../.. 
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Proof of Proposition 10.1. Using Lemmas 10.4 and 10.5 (with k = 1), we see tliat in 
o"„2 one may replace the term em(-^)"~^ with e~* at the expense of an error, not exceeding 

n-eJl + \tD • (l + \t\^^') e'^'/^ < <1_ h + |i|4m2) ^-'V^^ (10.7) 

where e'^^ 0. The same is true for the first m derivatives of cr„2- 

Now consider the products yk{t) = [{v{^) - em{^))'^ em{^)"' ^ appearing in (7„3 with 
2 < A; < n. Writing 

and combining Lemmas 10.3 and 10.5 (which give estimates that are independent of j), we get 
\yi\t)\ < 2f •£„(l + |t|^"')n-^'=/2.C(l + |t|2"^')e-("-'=)*'/(2") 

< (i , UiWn g_(n-fc)t2/(2n) 

- ^sk/2 V I I y 

n 
k=2 

n 
k=2 

For s > 2, we have 6n = n^*/^ g*^/2n _ q^i^-^^ uniformly in the interval \t\ < cn}/^ . So, 

(1 + 5nT = e"^''g(i+^") = e"('5n+0('5^)) = i + nSn + ^ {n5n)^ + nO{dl). 
Hence, for all n large enough, the expression in the large bracket in (10.8) does not exceed 

lM„y^+o(«)<^+o(^). 

It remains to compare this bound with (10.7) and (10.3), and then we arrive at (10.2). 

Finally, in the case s = 2, the expression in the large brackets in (10.8) is uniformly bounded 
in \t\ < cn}l^ . Thus Proposition 10.1 is proved. 



Therefore, 



11. Liouville Fractional Integrals and Derivatives 

In this section we recall basic definitions and some results on Liouville fractional integrals 
and derivatives, and refer to [S-K-M], [K-S-T] for proofs and a more detailed exposition. At 
the end of the section we also formulate some special estimates for such operators. The proof 
of Proposition 11.3 below is rather routine and is therefore postponed to the next section. 
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Let a denote a real number with < a < 1, and let y = y{t) denote a (measurable) 
function defined for t > 0. The Liouville left- and right-sided fractional integrals on the 
positive half-axis R"*" = (0, -|-oo) of order a are defined by 



it - xf-" 



{x > 0). 



The equalities are understood in the usual way (as Lebesgues integrals), if y is sufficiently 
"nice". According to a theorem by Hardy abd Littlewood, /q_|_ and are extended and act as 
bounded linear operators from LP(R+) to L'^(R+), where 1 < p, (/ < +oo, if and only if p < ^ 
and q = ■ They represent particular cases of the so-called Liouville (or Riemann-Liouville) 
fractional calculus operators. 

The Liouville left- and right-sided fractional derivatives on the positive half-axis are defined 

by 

{D^+y){x) = ^ (/i;"j/)(x), {D-y){x) = ^ (/i-"y)(x) {x > 0). 

The equalities are valid for sufficiently "nice" functions, including the class C^(Ti'^) of all 
infintely differentiable functions on R"*" with a compact support (which can be used to approx- 
imate functions from larger spaces). 

For example, for any complex number A, such that Re(A) > 0, 

(/^e-^*)(x) = e-^^, {D^e-^'){x) = A°^ e"^^, (ILl) 

where the principal value of the power functions is used. 

We cite two standard facts about these operators (see [K-S-T], p. 75 and p. 83). 

Proposition 11.1. For all sufficiently "good" functions y on R"*", 

iD^+I^_^y)ix) = yix), {D^I^y^x) = y(x). 

The equalities are extended to the space L^(R+). Moreover, if additionally y{x) = o(x") for 
X 0, then 

{lS+D^+y){x) = y{x). 

Define the linear spaces Iq_^_{LP(R^)) and /"(L^'(R"'")) as the images of LP(R'^) under the 
operators and 7", respectively. 

Proposition 11.2. For all sufficiently "good" functions f and g on R"'", 

f +0O 

g{x){D^_f){x)dx. (1L2) 



r+oo /-H 

/ f{x){D^^g){x)dx= / 
JO Jo 



The equality may be extended to all f G I^{U'(R'^)) and g G Iq^{L'^(R'^)) withp,q > 1, such 
that i i = 1 -h a. 

This is a formula for fractional integration by parts. 
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Now, let V be a function of bounded variation on the real line, also viewed as a finite 
measure, and denote by |F| its variation (as a measure). Define the Fourier-Stieltjes transform 

/ + 00 
e^*^ dV{t) X G R. 
-oo 

For our purposes, the following proposition will play a crucial role in the study of the local 
limit theorem with fractional moments. 

Proposition 11.3. Let g{x) = V(x) h{x), where h{x) is a continuously differentiable 
function on the real line, such that, for a given integer m > and < a < 1, as \x\ —^oo, 

\h{x)\ + \h'{x)\ = 0(1x1-^^+"*+")). 

// 1^1"^+" d\V\{t) < +00 and F(^)(0) = 0, for aUk = 0,...,m, then {Dg^g){x) exists for 
all X > and satisfies with some constant C , independent of x, 

\{D-^g){x)\ < J ^ min{\u\,\ur'}\ur\V\{dn). (11.3) 

Here a' = a in case m = 0, and a' = in case m > 1. In addition, for all t real, 

/ e'^'^ {D^+g){x) dx = i-it)" e'^^g{x)dx. (11.4) 
Jo Jo 

More precisely, the Gaussian function h{x) = e~^^/^ and its derivatives will only be needed 
in this proposition to obtain the desired decay for the inverse Fourier-Stieltjes transform. 



Proposition 11.4. For all functions V and h as in Proposition 11.3, for all t real, 

+00 

e**^ V{zx) h{x) dx 
where £{z) is hounded in \z\ < 1 and satisfies e{z) — >■ 0, as z — > 0. 



< ' I 1^ e(z), (11.5) 
- (l + |i|)" 



Proof. Let < 2; < 1. We apply Proposition 11.3 with the function Vz{u) = V{u/z) 
in place of V in which case Vz{x) = V{zx). Then, for the function gz{x) = V{zx)h{x), the 
fractional derivative {DQ_^gz){x) exists for all x > and satisfies (11.3). 

In order to unite both cases, use l^l" < 1^1" for \u\ > 1 and write (11.3) in a slightly weaker 
form 

m+9.)ix)\<j^^Siz), (11.6) 

where 

/+00 
min{|zn|,|znr} \ur\V\{du). 
-00 

This integral is finite and behaves like o{z"), as z — )■ 0. Indeed, split it into the two integrals 
in terms of the (finite positive) measure W{du) = \u\"^ \ V\{du) as 

d{z) = z'^Ioiz) + z'^hiz) = / dW{u) + / dW{u). (11.7) 

J\u\<l/z J\u\>l/z 
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By the moment assumption on V, we have j dW{u) < +00, so Ii{z) — )• 0, as z — )■ 0. As 
for the first integral, note that z^^^l^l < \u\°' in the region \u\ < 1/z. Hence, the functions 



fz{u) 



A-al 



1{|m|<i/2} have an integrable majorant f{u) = with respect to W. Since 



also fz{u) — 0, as z ^ 0, one may apply the Lebesgue dominated convergence theorem, which 
gives Io{z) = f fzdW ^0. Thus, from (11.6)-(11.7), 



m^9z){x)\<j^^e{z), 



(11.8) 



where e{z) = — )• 0, as 2; ^ 0, and supo<2<;^ e{z) < +00. 

Now, using the bound (11.8) in (11.4), we get | ff^°° e^'^^ gz{x) dx\ < ^e{z). Obviously, a 
similar inequality will hold as well when integrating over the negative half-axis. Therefore, 

f+00 ^ 

e'*''gz{x)dx 



f 

J — ( 



< r^e(z). 



This estimate implies (11.5) in case of large values of say, when \t\ > 1. The remaining 
range |i| < 1 can be treated by straightforward arguments. 

By the assumption on the decay of h, its Fourier transform h{t) = e**^ h{x) dx is 
well-defined, bounded, and has bounded continuous derivatives up to order m + 1. Introduce 
Taylor's approximation for h up to order m at a given point t, i.e., the function 



{Smh){t,u) =2^ ^, 
k=0 



U 



u eR. 



Prom Taylor's theorem it follows that 

\h{t + u)- {Smh)it,u)\ < Mmin{|nr,|'ur+^} 
with some constant M independent of t. Now write 



(11.9) 



/+00 r+00 
e^^"" gz{x) dx = I e'^''V{zx)h{x)dx 
-00 J —t 



' —00 



+00 



h{t + zu) dV{u). 



The assumption V{(}) = ■■■ = y("*)(0) = implies that /+^(S„/i)(t, u) dV{u) = 0, for all t. 
Therefore, using (11.9), we finally get 



\9zm 



{h{t + zu) - {Smh){t, zu)) dV{u] 



/+00 
mui{\zur ,\zur+^} d\V\{u) = o{z'^+'^). 
-00 



Note that the last relation has been already discussed in the previous step. 
Thus Proposition 11.4 is proved. 



Remark 11.5. The second part of the above proof also covers the limit case a = of the 
inequality (11.5), which may be written as 

f+00 

(^^''V{zx)h{x)dx < \zY^e{z). (11.10) 



/" 

J — ( 
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More precisely, for this assertion we only need the assumptions {ul™" \V\(du) < +00, 
V{0) = ... = fM(o) = 0, and 1x1"^+^ dx < +00. Proposition 11.3 is irrelevant in 
this case. 



12. Fourier Transforms and Fractional Derivatives 



In this section we give the proof of Proposition 11.3. By its very definition, 



r(a) {D^+g){x) 



d 



hit) 



V{t) dt, 



dx Jo (x - t)°' 

provided that the derivative exists (where < a < 1). 

Given m > integer, introduce the function of the real variable 



(12.1) 



k=0 



The assumptions on V imply that the first m moments of the measure V are vanishing, i.e., 
X!^ dV{u) = 0, for A; = 0, . . . , m. Hence, 

/oo 
r^m{tu)dV{u). 
-00 

Respectively, changing the variable in (12.1) and applying Fubini's theorem, one may write 



/+0O r i-x 
/ h{x - t) rim{{x - t)u) 
-00 L-'O 



dt 



dV{u). 



d^ /•+00 I _ ^ ^ ^ dt 

dx J-00 VJo ' ' ' ' 

Wc intend to move the differentiation inside the outer integral. To justify this step, consider 
the derivative with respect to the inner integral, 

rx 



-^(3;,?^) = ^^ h{x-t)r]m{{x-t)u)^. 



Lemma 12.1. Let h{x) denote a continuously differentiable function on the real line, such 
that \h{x)\ + \h'{x)\ = 0(a;-(2+"^+")), as \x\ — )• 00. Then, for all u and x > 0, 

\Iix,u)\ < C(l + x)-(^+") min{|^x|, (m = 0), 

\Iix,u)\ < C(l + ar)-(^+") min{|tiP,|-ur+i} (m > 1) 

with some constant C, depending on h and a, only. 



Proof. Put ^u{t) = h{t)r)miut), so that 

I(x,u) = f £„(x — t)— . 
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In this case we may interchange differentiation and integration. Thus, using -g^^uix — t) = 
— ■§liu{x — t) together with //^(O) = 0, we can write 

Iix,u) = -J^ C'ui^-t)-. 
Assume that x > 1 and spht the integration into the two regions such that 

I{x,u) = Io{x,u) + h{x,u) = - i'u{x-t)^- C{x-t)^. 

The integral Ii. 

Integrating by parts, we have 

h{x,u) = aj^ - *) ^ - - !)• (12-2) 

To analize this integral, we use the elementary bound 

|?7„(t)| <4min{|C,IC^'}, t € R. (12.3) 

Indeed, from Taylor's formula it follows that |?7m(i)| < (m+iy. ' "^^^^ settles (12.3) in case 
\t\ < 1. In the other case \t\ > 1, just write 



m I |fc ^ ™ 1 \ 

\^m{t)\ < 1 + E IT ^ 1*1"^ 1 + E ^ < (1 + ^) 1*1 
fe=0 ■ ^ k=0 



thus proving (12.3). 

This bound implies that 

\iu{?^-t)\ < 4|/i(x-i)|min{|(x-i)ur,|(x-t)ur+^}. (12.4) 

By the assumption on h, we have \h{x — 1)| < Cx~*^^+"^+"), so, by (12.4), 

\iu{x-\)\ < Cx-(^+")min{|wr,|nr+i} (12.5) 

with some constant C. 

In the region 1 < t < x\ = we use the bound \h(x-t)\ < Cx-^^+^+a) ^i^^ a constant 
independent of t and x. Hence, by (12.4), in this region 

Mx-t)] < Cx-(^+°^Vin{|nr,|-ur+^}, 

and 

\Ux -t)\-^< Cx-(i+-) min {|nr, \ur+^}. 
In the second region xi < t < x, just use < Cx-(i+"). Then, by (12.4), 
/ \Ux-t)\-^ < 4Cx-(i+")H-/ \h{x-t)\{x-trdt 

J X\ J Xl 

f{x-l)/2 

Jo 
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since the last integral is uniformly bounded. Similarly, again by (12.4), 

Collecting the bounds for the two regions, we get 

\Ux-t)\ < Cx-(i+")min{|nr,|t.r+i} 

with some constant C. Applying it together with (12.5) in (12.2), we arrive at 

h{x,u) < Ca;-(i+")min{|ur,|«r+i}. (12.6) 

The integral Iq. 

Now, let us turn to the integral Io{x, u) = — ^'u{x — t)^. After differentiation and using 
the identity r/^^ = irjm-i (with the convention that r)-i{t) = e**), one may represent it as 
-^0 = Iq,i +Io,2, where 

(It 

Iq,i{x,u) = -iu J r]m-iiix-t)u)h{x-t)—, 
7o,2(a;,n) = - rimiix - t)u) h' {x - t) — . 

Jo 5 

By (12.3), since x — t < x, 

|/o,2(a;,u)| <4 / \h'{x-t)\ mm{\xur,\xur+'}-. 

Jo * 

Using the assumption h'{x) = 0(x~'^^+'"+"^)) , we get 

\Io,2i^,u)\ < Ca;-(^+") mm{\ur,\ur+^} (12.7) 

with some constant C. 

As for the integral /o,i) first rewrite it as Iq^i{x,u) = —iu Jq h{x — t)d(jn-i{t), where 

/- r+\ rim-iiix - w)u) 
(m-iit)= / dw. 



Integrating by parts, one may represent it as /o,i = -^0,1,1 + -^0,1,2, where 

Io,i,i{x,u) = -iuh{x -l)Cm-i{^), 



Io,i,2i.x,u) = -iu h'{x -t)Cm-i{t)dt. 
Jo 



Claim. For all X >l,t& [0,1], we have \(^-i(t)\ < Cmin{l,|tt|" ^} with some constant C, 
while in case m > 1, 

\Cm-i{t)\<Cx^min{\ur-\ \ur}. 



Proof. If m = 0, whenever u ^ 0, 

dw 



IC-iWI 



Jo u;" ' ' Jo 



<C\u\'^-^. 
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On the other hand, \C-i{t)\ < /J ^, so \C-i{t)\ < Cmin{l, 

If m > 1, we just appeal to the estimate (12.4), which, for all w G (0, 1) and x > 1, yields 
\r]m-i{{x - w)u)\ < Amin {{x - w)""-^ \ur-\ {x-w)"'\u\"'} 
< 4x™min{|ur-M«r}. 
It immediately implies the desired estimate. 

Continuation of the proof of Lemma 12.1. Now, by the assumption on h, and using 
the claim in case m = 0, we obtain 

\Io,i{x,u)\ < \Io,i,i{x,u)\ + \Io,i,2{x,u)\ < Cx-(^+")min{|u|,|nr}. 

Similarly, in case m > 1, 

\Io,i,i{x,u)\ + \Io,i,2{x,u)\ < C\u\ {\h{x - 1)\ + \h' {x-l)\)-x"' mm {\ur-\ {uD 

< Ca;-(^+°^Vin{|ur, |nr+^}. 

Thus, 

\Io,i{x^'^)\ < Cx-(^+"Vin{|M|,|M|"} (m = 0), 

|/o,i(x,tx)| < Cx-^^+"'^mm{\ur,\ur+^} (m > 1). 

Taking into account (12.7) and (12.6), we arrive at similar bounds for Io{x,u) and I{x,u), 
which arc equivalent forms of the desired bounds in the lemma in case x > 1. 

Finally, let us only note that the case < x < 1 may be treated in a similar manner (with 
simpler estimates). Thus Lemma 12.1 is proved. 

Proof of Proposition 11.3. Finally, we are prepared to justify the differentiation step. 
Define 



dV{u), 



^{x)= / / Ux-t)- 

J-oo UO ^ . 

where, as before, ^u(t) = h{t)r]m{tu). Given x > and — >■ (with £„ 7^ 0, x + £„ > 0), write 

^(X + En) - i^ix) '■+°° r 1 rx+e. 



— / I{y,u)dy 



dV{u). 



C minllul, \u\°'} 



< C min{|u|, \u\°'}. 



J —00 \_ IT' <^ X 

By Lemma 12.1 in case m = 0, the expression in the square brackets is bounded in absolute 
value by 

1 fX+en 

- (l + y)-(i+")dy 

J X 

On the right-hand side the function is integrable with respect to the measure \V\, according 
to the moment condition on the function V . A similar conclusion holds in case m > 1 (with 
appropriate modifications in the estimate). Therefore, one may apply the Lebesgue dominated 
convergence theorem, which gives 

ip yx) = lim I I[x,u)aV{u). 
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Thus, the fractional derivative 

PoV)(^) = = I{x,u)dV{u) 

is well-defined for all x > 0. Moreover, from Lemma 12.1 we also obtain that 

/ + 00 
uAn{\u\^\\u\'^^}d\V\{u), 
-oo 

where m\ = 1, ra2 = q in case m = 0, and m\ = m, m2 = m + 1 in case m > 1. 

This proves the first assertion and the inequality (11-3) in Proposition 11.3. For the second 
assertion, apply Proposition 11.1 and the formula (11.2) for the fractional integration by parts 
with the functions /(x) = e~(^~'*)^ {e > 0) and use the second formula in (11.1) with X = e—it 
for the fractional derivatives of /. Then we obtain 

r+oo f+oo 

/ e-(^-^*) ^ (D^+g) (x) dx = {e- it)'' / e'^^''^^ ^ g{x) dx. 
Jo Jo 

Letting e — t- and using the intcgrability of both g and D^^g (due to (11.3)), we arrive in the 
limit at the required equality (11.4). Thus Proposition 11.3 is proved. 



13. Binomial Decomposition of Convolutions 

We shall now the probability densities p„ in Theorem 1.2. The next procedure is known; 
a related approach has been used, e.g., in [S-M], [I-L] to study the central limit theorem with 
respect to the total variation distance. 

Let < c < 1 be a prescribed number, m = [s], and n > m + 2. 

Without loss of generality, one may assume that no = 1, that is, = Xi has a density, say, 
p which may or may not be bounded. For definiteness, assume it is (essentially) unbounded, 
so that the integral 

b = p{x) dx 

J{p{x)>M} 

is positive for all M > 0. We choose M to be sufficiently large to satisfy, e.g., < 6 < | which 
implies 2n"^~^^ yn-m-i ^ ^n^ £qj, g^^j n > ni large enough. 
Consider the decomposition 

p{x) = ap{x) + bq{x), 

where a = 1 — b, and p{x), q{x) are the normalized restrictions of p to the sets {p{x) < M} 
and {p{x) > M}, respectively. Hence, for the convolutions we have a binomial decomposition 

n 
k=0 

Then split the above sum into the two parts to get p*^{x) = Pn{x) + qn{x), where 

n m+1 

Pn= Yl C^oJ'b'^'^ p*^ * qn=Yl C^a^b'^-^p*^ * 

k=m+2 k=0 
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Note that 



r+oo 

qn{x) dx = 

k=0 

f +00 



/+0O m+1 „ 

qn{x) dx=^C^ a^""-^ < n"^+^ bn-m-l ^ > 
-00 2 

Finally define 



■00 



/n 

Let us recall that Pn{x) = ^Jn p*'^{x^/n^ has the characteristic function Vnii) = ^("^)") 
where v is the characteristic function of X\. By construction, the densities p„ are bounded 
and provide a strong approximation for p^. Namely, we immediately obtain: 

Lemma 13.1. For all n > ni, \Pnix) — Pn{x)\ dx < c". In particular, for all t G R, 

\Vn{t)-Vn{t)\<c''. 

A similar inequality also holds for the first m derivatives of Vn and Vn with n large enough. 

The last assertion of the lemma needs a more detailed explanation which we postpone to 
the end of the section. 

We will also need some integrability properties for Vn and their first m, derivatives that are 
due to the boundedness of the probability density p{x) and the finiteness of the m-th absolute 
moment of X\ . 

Lemma 13.2. Let E|Xi|™ < +00, m > 2. There exist positive constants A and a, 

depending on Xi and m, such that, for all < T < \fn. 



/ 



\vn(t)\dt < Ae-''''^\ (13.1) 

/{|t|>T} 

A similar bound is also true for the first m derivatives ofvn {with arbitrary n> m + 2). 
Proof. Let p, q denote the Fourier transforms of p and q, respectively. Then 

By the Riemann-Lebesgue theorem, 

sup|p(t)|<7, sup|g(t)|<7 (0<7<1). (13.3) 

\t\>i \t\>i 

Hence, from (13.2), for all |t| > ^/n, 

I^"^*^'<l^l^(;^)l' ^ C^-'b--'<r-'\p{^)f. (13.4) 

In addition, by the Plancherel theorem, and using p{x) < M/a, we have 

\p(-=)fdt = 2TrVE p{xfdt<-—^. (13.5) 
-00 V" i-00 « 
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Therefore, integrating the inequahty (13.4), we get 



27rM 



/ \Vnit)\dt< 



■7 



,n-2 



^/n. 



(13.6) 



On the other hand, since both p and q represent the densities of probability distributions 
with finite second moments, their characteristic functions near zero satisfy 

|^(t)|<e--'*', \m\<e-''"'" (|i|<l) (13.7) 
with some constant a > 0. Hence, for \t\ < ^/n, (13.2) gives the estimate |S^„(f)| < e^'^^*^ and 



[ \Vnit)\dt < [ 

J iT<\t\<./n\ JT 



* dt < -e""^ . 



'{T<\t\<V^} 

Together with (13.6) the latter yields 



T<|t|<VS 



m>T} 



\vn{t)\ dt<- e-'^'^' + 7"~'V^- 



27rM 



Finally, since for the values < T < ^Jn one always has 7"^ < ^1 e '^^'^^ with some 

constants Ai and ai > 0, the desired bound (13.1) easily follows. 

As for the derivatives, a bound of this type can be proved by similar arguments, so let us 
restrict ourselves to the basic case of the m-th derivative (needed for the proof of Theorem 1.2). 

The condition E|Xi|™ = J \x\"^ p{x) dx < +00 implies a similar property for the densities 
p{x) and q{x). Hence, p{t) and q{t) have continuous derivatives up to order m, bounded in 
absolute value by some common constant. 

In view of (13.2), vli^\t) represents a linear combination of the terms 



d" 



dV 



^ I t \k^/ t \ n—k 

\/n' \/n' 



""dr 

r=0 



df"' 



^1 t \n—k 



(13.8) 



with integers m + 2 < k < n. Given < r < m, by the chain rule (2.3), the r-th derivative of 
the function p(-^)'= represents a linear combination of the terms 

\f\..p(-){-Lf^ (13.9) 



n 



-r/2 -f ^ \k-{ki-\ hkr)^// 



P 



n 



n 



over all integer tuples (ki, . . . , kj.), such that ki + 2k2 + • • • + rkr = r and ki + ■ ■ ■ + kr < k 
(kj > 0). Moreover, the coefficients in that linear combination do not depend on n, and the 
total number of such terms is bounded by a quantity that depends on m, only. 

Using ki + ■ ■ ■ + kr < r < m and the boundedness of the derivatives, the absolute value of 



the expression (13.9), as well as the sum of all such terms is bounded by |p(;^j 
constant factor. 



\k—m 



up to a 



Similarly, the (m— r)-th derivative of the function ^(-^)" represents a linear combination 



of the terms 



n 



(m-r)/2 ^ ^(n-fc)-(fciH hkm-r) ( J_\ki 



q{m-r) ij_\kr, 



(13.10) 



over all integer tuples (fci, . . . , k^-r), such that ki + 2^2 + • • • + (m — r)km-r = m — r and 
ki + ■ ■ ■ + km~r < n — k {kj > 0). Again, the coefficients in the linear combination do not 
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depend on n, and the total number of such terms is bounded by a quantity depending on m, 
only. Since ki + ■ ■ ■ + km-r < min(n — k,m — r) < min(n — k,m), the absolute value of the 
expression (13.10), and the sum of all such terms is bounded by |^(_^)|("'-*;)-min(n-fc,m) 

a constant factor. 

Thus, (13.8) is bounded in absolute value by 

with some constant C, depending on Xi and m, only. It then follows from (13.2) that 

C f-ik „kun-k \^( * ^\k-m \ t | (n-fc)-min(n-fc,m) 

Now, like in the previous step, using (13.3), for all \t\ > \fn^ we get 

W'^ /-(fe fetn-fe (n-m-2)-min(n-fe,m) 



< C_ ^ c^'a*^6-*^|p( * )|'=-"^ |g(^)|('^-'=)-'"^'^("-'=''"). (13.12) 



fc=m+2 

Integrating this inequality with the help of (13.5), we obtain that 

/ |4-)(t)| dt < 7"-2— ^v^. (13.13) 

J{\t\>^} « 

In addition, using (13.7), for |t| < ^/n, the product (13.11) is bounded by Ce"'^'^^*^, where 

d = — ((ra — m) — min(n — A;, m)) > d' = — ((n — m) — min(n — m — 2, m)) . 

If n > 2m + 2, then d' = ^^^^^ > In the other case m + 2 < n < 2m + 2, we also have 

d' = ^> Hence, (13.11) is bounded by (7e-^'*'/("*+^), and we derive from (13.12) 

C 



|4'»)(t)| < C^a'=6"-'=e-'^ * /('"+^) < Ce-'^ * 

fc=m+2 

It remains to integrate this inequality to get 

/ \vt\t)\ dt<C [ e-'^^t- lira^X) < ,-a^TV(m+l). 

J{T<|t|<Vn} JT<\t\<s/n ^ 

Together with (13.13), it yields the desired estimate (13.1). Thus Lemma 13.2 is proved. 

Remark 13.3. If the density p is bounded, the decomposition procedure is not needed, 
and then Lemma 13.2 should read as follows. Let E < +00, m > 2, for a random variable 
having a bounded density. There exist constants A and a > 0, such that for all n > 2, 



t 



\v{-^)rdt<Ae-'''^ , 0<T<V^, (13.14) 
{|t|>r}' 'V^ ' 
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where v is the characteristic function of X. A similar bound holds as well for the first m 

t \n 



derivatives of v{^)'" with arbitrary n> m + 2. 



Proof of Lemma 13.1. By the construction, for all n > rai, 

f +00 



/ 

J — ( 



\Pn{x) - Pn{x)\ dx < 2/3n < c"*, 



SO \vn{t) — Vn{t)\ < c", as well. In order to extend this inequality to the derivatives, recall the 
representation (13.2) to write 

Vn{t) - Mt) = ^1 - (13.15) 



where 



So = 



m+1 



As before, we will only consider the case of the m-th derivative. 

It has been shown in the proof of Lemma 13.2 that, given m + 2 < k < m, the function 
pi^^Y ^{^)^~'^ the m-th derivative, bounded in absolute value by the expression (13.11). 
So, it is bounded by a constant C, depending on Xi and m, only. In the general case, including 
the values < < m + 1, (13.11) should be replaced with 

^1^1- t , imax(A;— m.O) i ^/ t , i (n— fe)— min(n— fe.m) 

which is also bounded by C. Therefore, from (13.15), 



4-)(t)-^(-)(t)| < Cj^aH^-' + CY'^n^'b--" 



" k=m+2 k=0 



= 2C/3„ < c", 

where the last inequality holds true for all n starting with a certain ni. 
Thus, Lemma 13.1 is proved. 



14. Proof of Theorems 1.1 and 1.2 

We are prepared to make the last step in the proof of Theorems 1.1 and 1.2. Recall that 

s > 2, m = [s], and put a = s — m. 

Let v{t) be the characteristic function of Xi and Vn{t) = t;(--^)" be the characteristic 

function of Sn- We will assume that all Sn have densities pn (since only minor modifications 
have to be done in the more general case, where Sn have densities for all n large enough). 
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If priQ and therefore all pn with n > no are (essentially) bounded for some no, then there 
is no need to use the binomial decomposition of the previous section, and we put pn = Pn- 
This case corresponds to Theorem 1.1. Otherwise, if p„ are unbounded for all n > 1, then the 
binomial decomposition is applied to p = pi, and we obtain the modified densities pn together 
with the associated characteristic functions which we considered in the previous section. 
Thus, the requirement c) in Theorem 1.2 is met. 

The inversion formula. 

The characteristic functions Vn have continuous, bounded derivatives up to order m, that 
are integrable according to the inequality (13.1) of Lemma 13.2, or (13.14) of Remark 13.3. 
Hence, by the inversion formula, 

1 r+QO 

{ixypn{x) = —j_ e-''^vlP\t)dt, p = 0,l,...,m. 

By the construction, the approximating functions ipm{x) = ^{x)+YlT=i Qk{x) n~^/'^ , which 
appear in the relation (1.3), have the integrable Fourier transform 

rrt-2 



( Hi — ^ \ 
k=l ' 



Consequently, for all p = 0, 1, . . . , m, 

{ixf {pn{x) - ^^(x)) = ^ y ^ e-^*-(iy(^')(t) - 4^}{t)) dt. (14.1) 

Our task is thus to give proper upper bounds on the absolute value of these integrals in the 
particular cases p = and p = m. 

What is rather standard, one should split the integration into two regions. Given T„ +oo, 
< r„ < -y/n (to be specified later on), let 



fn,p = / e-*- {v^^ {t) - u^) it)) dt, Jn,, = f e-*- {v(P^ it) - it)) dt. 

J\t\<Tr,. J\t\>Tr, 

\t)\ dt < Ae~ 



/ 

J\t\ 



l\t\<T„ J\t\>Tn 

It should be clear that 

'\t\>T„ 

with some positive constants A and a, depending on m. By Lemma 13.2 and Remark 13.3, 
we have a similar bound for Vn\t), whenever n > m + 2, so 

kn,p| <^e-<^'^". (14.2) 

The integral 

To treat this integral, we subtract and add Vn\t) inside the integrand and apply Lemma 
13.1 (the second part). Then it gives (using T„ < ^/n) 

\In,p\ < \In,p\ + c"^/n, 

where < c < 1 is the prescribed parameter in Theorem 1.2, and 



In,p= f e-^'^{v(f{t)-u!^\t))dt. 
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Using (14.1)-(14.2), we obtain that, for all x, 



up some positive constants A and a. 



(14.3) 



Proof of (1.4) in case |x| > 1. 

Note that for |a;| < 1, the relation (1.4) follows from (1.3). As for the values \x\ > 1, only 

the value p = m is of interest in (14.3). 

The integral I'^p can be treated with the help of Theorem 1.3. It gives that in the interval 
\t\ < ci'n}/^ with some constant < ci < 1 we have 



where the remainder satisfies 



|r„,| < e" 



C 



+ 



n 



(m-l)/2 ^5-2 / • 



Here C, c and e„ are some positive constants, such that e„ — t- 0, as n — t- oo. Hence, assuming 
that Tn < cin^/^ and noting that cP'^fn will be absorbed by other remainder terms, we get 
that 



where 



c 



+ 



I" 



-itx 



\t\<Tr, 



dtp 



n{m-l)/2 nS-2' 



n 



dt. 



(14.4) 



(14.5) 



Now, one can differentiate inside the last integral, which will lead to the terms, containing 



up to polynomial factors (due to the property that v has m bounded derivatives). Hence, 



integration in (14.5) may be extended to the whole real line at the expense of an error not 
exceeding Ce~^"/^. Hence, (14.4) may be replaced with 

C 



'■n,p 



+ 00 



-itx 



27r 

(F_ 

dtp 



+ 



(14.6) 



-t^/2 



dt. 



„(m-l)/2 „s-2' 

where 

Letting p = m and w{t) = v{t) — em{t), and performing differentiation, rewrite the above 
integral as 

In,m = E # / -^'^ (-^) hm-kit) dt, (14.7) 

fe=o " ' -^-^ 

where h^-^(t) = (e-*V2){m-fe) ^ (_i)--fc^^_^(i) e-tV2. 

Recall that w{t) = V{t) represents the Fourier transform of a finite signed measure, V, such 



that /_ 



|F|(ciit) < +00 (where \V\ denotes the variation of V, treated as a positive 



finite measure). In addition, the first m derivatives of w are vanishing (cf. Section 7 and 
Proposition 7.4). Hence, w^''\t) = Vk{t) represents the Fourier transform of a finite signed 
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measure 14, such that /_!^ |^(|(™-^)+" |Vfc|(d?x) < +00, and also the first m — k derivatives 
of Wk are vanishing. Therefore, we are in position to apply Proposition 11.4 to the functions 
Vk{t) in place of V, h^-k in place of /i, and with m — A; in place of the parameter m. Choosing 
z = l/-\Ai, the inequalities (11.5) and (11.10) (cf. Remark 11.5) give 



/_ 



+00 J. 



n' 



< ^ fl + bir" 



with some sequence £„ — > 0, as n — >■ 00. Applying this bound in (14.7), we obtain 



and then (14.6) with p = m yields 

|xr |p.(x) - ^.(-)| < + \xr{Ae-'^'^^ + + . (14.8) 

It remains to involve information about the possible growth of T„. But, as we have seen, 
one could choose r„ to be of order n~^/^ regardless of s. With this choice (14.8) leads to the 
announced inequality (1.4) of Theorem 1.2. 

Proof of (1.3). 

Let us return to (14.3). The integral I'^p can also be estimated by virtue of Proposition 
9.1 and Propositions 5.1-5.2 (cf. Corollary (9.2)). Namely, they give that 

t;(f)(t)-w(P)(t) = o(n-('*-2)/2)e-*'/^ p = 0,l...,m, (14.9) 

uniformly over all t in the intervals \t\ < Tn, where T„ are of order n~^^^ in case s > 3 and 
of order n^i^^'^yC^^) in case 2 < s < 3. If s = 2, we may only have T„ +00. Clearly, in all 
cases I'j^p = o(ra~(*~^)/^), and (14.3) yields 

\x\P \pn{x) - <fm{x)\ < o(n-(^-2)/2) +^e-'^'^«. (14.10) 

It remains to apply this inequality with p = and p = m. 
Theorems 1.1-1.2 are thus proved. 

Remark. If E < +00, m > 2, but (pm are constructed with the help of the same 

cumulants 73, . . . ,7™, (like in the case m < ,s < m+ 1), the relation (1.3) for both Theorem 1.1 
and Theorem 1.2 may be sharpened. Indeed, by Proposition 9.1 (second part), (14.9) should 
be replaced with a stronger relation 

v^r?\t) - ut\t) = 0(n-(™-i)/2) e-*V^ p = 0, 1 . . . ,m, 

which holds uniformly in the intervals \t\ < civ}/^. Respectively, it provides a stronger variant 
of (14.10), namely, 

(1 + \pn{x) - Vm{x)\ = 0(n-(— 
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